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ABSTRACT 

A  statistical  survey  of  the  wind  veering  in  the  lowest  two  kilometers 
is  made  for  the  oceans  of  the  Northern  Hemisphere.  The  data  sample 
consists  of  over  100,000  pibal  and  rawin  soundings.  This  includes  ob¬ 
servations  made  from  ships  in  the  0^-70°  latitude  belt  and  observations 
at  atolls  and  small  islands  in  the  5° - 3 5CT  latitude  belt.  The  wind  veer- 
ing  is  stratified  by  latitude,  season,  wind  direction  and  wind  speed.  T+re 
results  show  average  veering  values  of  8-12  degrees  in  the  lowest  km 
layer  and  0-3  degrees  in  the  second  km  layer.  The  average  veering  shows 
little  latitude  or  seasonal  variation.  The  veering  increases  with  wind 
speed  in  the  middle  latitudes  but  not  in  the  tropics.  As  expected  from 
thermal  wind  considerations,  south  winds  show  a  characteristic  veering 
of  20°  more  than  north  winds.  This  thermal  wind  influence  is  largely 
eliminated  by  assuming  constancy  of  thermal  wind  with  height  and  sub¬ 
tracting  the  second  km  layer  veering  from  the  first. 

It  is  observed  that  the  height  of  the  planetary  boundary  layer  does 
not  increase  towards  the  equator  and  can  be  specified  by  vertical  stabil¬ 
ity  and  turbulent  intensity  considerations  without  resort  to  rotation  argu¬ 
ments.  The  frictional  induced  kinetic  energy  (KE)  dissipation  is  substan¬ 
tially  greater  than  the  cross  isobaric  KE  generation.  In  the  tropics  the 
dissipation  to  generation  ratio  is  more  than  five  to  one.  The  oceanic 
boundary  layer  can  be  maintained  only  by  a  downward  transport  of  KE 
from  higher  levels. 


I.  BACKGROUND 


Few  meteorologists  doubt  the  crucial  importance  of  the  planetary 
boundary  layer  (p.  b.  1. )  in  atmospheric  dynamics.  Indeed,  a  large  por¬ 
tion  of  the  Stockholm  GARP  Report  (1067)  was  devoted  to  discussions  of 
the  planetary  boundary  layer  (p.b.l. )  by  Charnock  and  Ellison,  Black- 
adar.  Monin  and  Zilitinkenich,  and  Priestly.  Sheppard  (1969)  has  given 
a  summary  of  our  present  p.b.l.  knowledge  with  regard  to  large  scale 
global  modeling. 

The  associations  of  cumulus  convection  with  positive  p.b.l.  relative 
vorticity  is  very  evident  along  fronts  and  on  the  cyclonic  shearing  side 
of  the  trade  winds.  This  cumulus  convection  is  thought  to  be  primarily 
induced  by  frictional  wind  veering  in  the  p.b.l.  Where  tropospheric  ver¬ 
tical  wind  shears  are  small,  and  other  conditions  are  favorable  non-lin¬ 
ear  feedbacks  and  flow  intensification  can  occur.  This  is  the  so-called 
CISK  mechanism  which  is  generally  believed  necessary  for  the  inten¬ 
sification  of  tropical  storms  and  cloud  clusters.  Over  the  oceans  the  top 
of  the  p.b.l.  (or  where  the  surface  induced  mechanical  wind  veering  stops) 
is  observed  to  be  close  to  the  base  of  the  cumulus  clouds.  These  oceanic 
cumuli  are  dependent  on  sub-cloud  layer  convergence,  a  significant  por¬ 
tion  of  which  is  believed  to  be  mechanically  forced  by  surface  friction. 

The  assumptions  of  the  Ekman  theory  (1905)  for  this  layer  (steady 

*Conditional  Instability  of  the  Second  Kind  (CISK).  Physical  idea  pro¬ 
posed  by  Ooyama  (1964)  and  Charney  and  Eliassen  (1964)  and  so  much  dis¬ 
cussed  in  connection  with  the  GARP  and  GATE  experiment. 
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motion,  no  advection,  neutral  stability,  no  thermal  wind  and  constant 
coefficient  of  viscosity)  are  almost  never  simultaneously  valid.  The 
question  is,  how  invalid  are  these  assumptions  in  the  real  atmosphere? 
How  much  does  the  so-called  "Ekman  Layer"  of  the  real  atmosphere  sub¬ 
scribe  to  Ekman's  theory?  Granted  that  Ekman's  theory  fails  for  the 
individual  sounding,  nevertheless,  is  the  concept  a  useful  one  in  a  sta¬ 
tistical  sense  where  the  thermal  wind,  stability,  advection,  and  tendency 
variations  are  mostly  averaged  out  of  the  data  sample? 

Turbulence  theory  is  not  developed  to  the  point  of  allowing  for  a  prac¬ 
tical,  objective,  generalized  p,  b.  1.  model  which  will  handle  all  of  the 
parameter  variations  of  the  individual  time  periods  in  terms  of  a  synop¬ 
tic-scale  data  sample.  Yet  future  global  numerical  forecasting  will  not 
be  very  successful  unless  a  realistic  treatment  of  the  p.  b.  1.  is  accom¬ 
plished.  Probably  we  will  not  be  able  to  wait  for  a  satisfactory  turbu¬ 
lence  theory  to  be  developed  before  we  proceed  in  our  attempts  to  para¬ 
meterize  the  effects  of  the  p.  b.  1.  in  terms  of  the  larger-scale  motions. 
For  these  reasons  the  author  feels  that  at  this  stage  our  knowledge  of 
the  p. b.  1.  will  be  most  rapidly  increased  by  going  directly  to  the  obser¬ 
vations  and  studying  empirically  how  all  of  the  measurable  parameters 
vary. 

Previous  Observational  Studies.  Most  boundary  layer  research  has 
concentrated  on  the  lowest  10  to  20  meters  of  the  atmosphere.  Thousands 
of  papers  have  been  published  on  this  lowest  atmospheric  layer.  In  a 
comparative  sense,  the  10  meter  to  1  km  or  Ekman  layer  has  been  largely 
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neglected- -especially  observationally. 

Previous  observational  studies  of  the  entire  p.b.  1.  frictional  veering 
layer  have,  in  general,  been  limited  in  location  and/or  length  of  record. 
Some  of  these  noteworthy  pioneering  observational  studies  of  the  p.  b.  1. 
include: 

Gregg  (1922)  -  4  Kite  locations  over  the  U.S.  (discussed  and  rework¬ 
ed  by  Johnson,  1962). 

Westwater  (1943)  -  North  Sea  and  Trade  Winds 
Frost  (1948)  -  A  number  of  locations  in  Europe. 

Gordon  (1950a,  1950b,  1952a,  1952b)  -  North  Atlantic 
Sheppard,  Charnock,  and  Francis  (1952)  -  Scilly  Islands 
Sheppard  and  Omar  (1952)  -  Trades 

Charnock,  Francis,  and  Sheppard  (1956)  -  Anegada,  Virgin  Islands 
The  latter  three  studies  made  use  of  multi -theodolite  observations 
to  obtain  higher  wind  accuracy.  This  often  allowed  an  estimate  of  wind 
component  correlations,  stress,  and  exchange  coefficients,  etc.  which 
contributed  to  a  much  better  layer  description. 

The  most  extensive  observational  work  on  the  p.  b.  1.  has  been  accom¬ 
plished  by  the  University  of  Wisconsin  group  led  by  H.  Lettau  [Lettau 
and  Davidson  (1957),  Lettau  (1959),  Lettau  et  al.  (1962),  Lettau  (1967), 

B.  Lettau  (1967),  plus  many  other  papers]  and  his  most  active  former 
colleague  E.  Kung  [1963,  1967,  1968,  and  othei  papers].  Other  recent 
very  worthy  observational  p.b.  1.  studies  have  been  accomplished  by 
Findlater  et  al.  (1966),  Mendenhall  (1967),  Bonner  and  Smith  (1967),  Clark 
(1970),  Janota  (1971),  and  Cattle  (1972). 
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Additional  recent  observational  studies  which  have  dealt  with  the  p.  b.  1. 
have  included  the  work  of  Atkinson  and  Saddler  (1971),  Ilasso  (1970), 

Hasse  and  Wagner  (1971),  Estoque  (1970),  Aagaard  (1969),  Angell  et  al. 
(1968),  and  Kraus  (1966).  These  are  about  all  the  researchers  who  have 
made  significant  detailed  observational  studies  through  the  entire  p.b.  1. 

Most  of  these  later  studies  have  been  directed  to  specific  p.  b.l.  phy¬ 
sical  questions  and  have  been,  of  necessity,  restrictive  in  areal  extent 
or  time  sample.  Some  of  these  observations  were  taken  over  land.  A 
general  deficiency  exists  with  regard  to  oceanic  observations. 

With  recent  availability  of  magnetic  data  tape  storage  at  the  National 
Climatic  Center,  Asheville,  N.C.,  many  of  these  observational  deficien¬ 
cies  over  the  oceans  can  now  be  alleviated. 

Use  of  Rawin  and  Pibal  Observations.  In  order  to  remove  data  sam¬ 
ple  restrictions  and  obtain  a  broader  view  of  the  p.  b.  1. ,  both  pil  al  and 
rawinsonde  data  are  used  in  the  present  study.  The  use  of  conventional 
pibal  and  rawinsonde  data  solves  the  problem  of  spatial  representation. 
These  observations  are  routinely  available  from  widely  separated  loca¬ 
tions  having  many  ditferent  environmental  characteristics.  The  more 
accurate  but  much  less  abundant  special  observations  are  not  used  (ex¬ 
cept  for  a  small  amount  of  1  967  Line  Island  Experiment  data  -  Madden, 
et  al.  1971).  A  sacrifice  of  individual  site  reprt  ^entativeness  for  a 
large  data  sample  felt  justified  by  the  fact  that  it  is  the  average,  broad- 
er-scale,  p.  b.l.  effects  which  are  being  sought.  A  large  data  sample  is 
more  likely  to  be  of  value  for  this  purpose  than  a  highly  accurate  data 


5 


sample  which  may  have  biases  characteristic  of  the  individual  site  ther¬ 
mal  wind,  stability,  etc.  Additional  observational  restrictions  include 
limited  vertical  resolution  and  real  turbulent  fluctuations.  The  author 
feels  that  these  latter  fluctuations  are  largely  random  and  are  eliminated 
in  a  large  data  sample. 

How  Much  p.  b.  1.  Wind  Veering  Is  There  over  the  Oceans?  Some  me¬ 
teorologists  say  that  the  Ekman  wind  veering  over  the  oceans  is  typically 
too  small  to  be  of  much  importance  and  that  consequently,  the  CISK  or 
low-level  frictional  forcing  idea  may  be  overemphasized  (Simpson,  1971). 
In  a  study  of  the  Trades  of  the  northeast  Pacific,  Riehl  et  al.  (1951)  ob¬ 
tained  a  mean  p.  b.l.  wind  veering  of  only  4°-6°.  What  is  the  typical 
frictional  wind  veering  over  the  oceans?  Many  other  fundamental  ques¬ 
tions  on  the  oceanic  p.  b.  1.  such  as  its  typical  height,  the  change  of  height 
as  one  approaches  the  equator,  the  typical  momentum  and  kinetic  energy 
dissipation  to  generation  ratios,  etc. ,  remain  to  be  answered. 

Let  us  first  compute  the  frictional  veering  through  the  lowest  few 
kilometers  at  many  global  locations  and  under  a  variety  of  conditions  be¬ 
fore  we  attempt  to  develop  a  generalized  theory  or  working  model  of  the 
p.  b.  1. 
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II.  OBSERVATIONAL  CC  NSIDERATIONS 

Veering  Variations.  The  turbulent  nature  of  the  p.  b.l.  limits  the  in¬ 
formation  that  can  be  gained  about  the  broader-scale  flow  from  individ¬ 
ual  observations.  Mechanical  eddies  on  gust  scales  of  50-500  m  can 
often  disrupt  the  upward  trajectory  of  the  balloon  and  typically  give  in¬ 
dividual  veering  values  unrepresentative  of  the  average  veering  condi¬ 
tions  surrounding  the  ascent  location.  Careful  double  or  triple  theodo¬ 
lite  observations  often  cannot  overcome  this  turbulence  induced  deficien¬ 
cy  of  individual  observations.  In  addition,  instrumental  inaccuracies 
occur  due  to  initial  balloon  tracking  pickup  from  ground  release  and  the 
long  time  averaging  interval  (1  to  2  minutes).  Extraneous  balloon  mo¬ 
tions  can  further  detract  from  the  representativeness  of  the  individual 
reports. 

Figs.  1  and  2  illustiate  the  large  6  and  12  hourly  1  km  veering  fluc¬ 
tuations  at  Weather  Ship  N  (80oN,  140°W)  and  at  Swan  Island  (17°N, 
84°W).  These  are  typical  of  ail  stations.  One  might  have  expected  a 
time-series  of  wind  veering  in  the  planetary  boundary  layer  tr  show  some 
degree  of  constancy,  especially  in  the  relatively  steady  trade-wind  re¬ 
gime.  This  is  not  the  case.  Observed  wind  veering  is  highly  variable. 
Figs.  3  and  4  show  a  typical  group  of  consecutive  rawin  profiles  of  ver¬ 
tical  variation  of  wind  direction  at  Wake  (20°N,  166°E)  and  Swan  Island 
(17°N,  84°W).  The  usual  autocorrelogram  of  lowest  km  wind  veering 
angle  for  increasing  time  lag  is  shown  in  Fig.  5.  At  6  hours  the  auto¬ 
correlation  is  C.3,  at  12  hours  and  beyond  it  is  zero.  Fig.  6 


Fig.  1. 
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Time-series  of  observed 
wind  veering  in  lowest  km 
layer  at  Ship  N  (30N,  ItOW) 
ob.ained  from  6 -hourly  ra- 
win  sov  ndings  (from  Men¬ 
denhall,  1967). 


Fig.  2. 

Time-series  of  observed 
wind  veering  in  lowest  km 
layer  at  Swan  Island  (17N, 
71W)  obtained  from  12- 
hourly  rawin  soundings. 
Dashed  lines  indicate  per¬ 
iods  of  missing  observa¬ 
tions  (from  Mendenhall, 
1967). 
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WIND  VEERING  RELATIVE  TO  2  Km.  HEIGHT 

Fig.  3.  Typical  individual  sounding  variation  of  rawin  measured 
wind  direction  in  lowest  2  km  layer.  Soundings  were 
taken  6  hours  apart.  Mean  veering  shown  by  the  solid 
line. 
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Fig.  6. 

Smoothed  frequency 
distribution  of  veer¬ 
ing  at  Johnston  Island 
in  the  layer  from  sur¬ 
face  to  900  mb.  Sample 
size:  3696  observations 
(from  Mendenhall,  1967). 
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UPPER  AND  LOWER  VEERING 
RELATIVE  TO  I  Km  LAYER 


Fig.  7.  Average  of  the  50  percent  largest  (or  upper  half)  and 
50  percent  smallest  (or  lower  half)  veering  values 
in  each  of  the  1st  and  2nd  km  layer.  One  km  wind 
is  used  as  reference,  heavy  line  is  the  mean  of  all 
data.  Veering  is  present  if  the  angle  shown  on  the 
abscissa  decrease  with  height. 
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shows  the  broad  frequency  distribution  of  1st  km  wind  veering  at  John¬ 
ston  Island.  This  is  the  normal  situation  at  all  stations.  This  veering 
is  due  to  both  turbulent  eddies  and  random  instrumental  errors. 

These  large  1st  km  veering  fluctuations  are  also  present  to  nearly 
the  same  degree  in  the  second  km  layer  as  seen  in  Figs.  7  and  8.  Here 
the  veering  data  was  ordered  from  highest  to  lowest  value  in  each  lay¬ 
er.  These  figures  show  the  average  of  the  50  percent  largest  (or  upper 
half)  and  50  percent  smallest  (or  lower  half)  veering  values  in  each  of 
the  1st  and  2nd  km  layers.  The  one  km  wind  is  used  as  reference.  Veer¬ 
ing  is  present  if  the  angle  shown  on  the  abscissa  decreases  with  height. 
Although  one  finds  a  similar  scatter  of  veering  in  both  layers,  the  av¬ 
erage  1st  km  veering  is  about  10°  more  than  the  2nd  km  veering.  Tur¬ 
bulent  component  corrections  must  be  higher  in  the  1st  km.  There  is 
no  indication  that  the  2nd  km  turbulent  and  instrumental  fluctuations  are 
not  random.  They  are  self  cancelling  in  a  long  period  average. 

Ship  Wind  Restrictions.  Wind  directions  are  least  accurate  at  very 
light  speeds.  Wind  speeds  are  least  accurate  at  very  strong  speeds. 

Ship  soundings  are  considered  less  reliable  with  elevation  angles  above 
60°  or  below  15°,  which  are  equivalent  to  average  wind  speeds  of  2-1/2 
and  19  m/sec  respectively  assuming  300  m/min  rate  of  rise.  Only  a 
small  percent  of  the  wind  speeds  were  below  2-1/2  or  above  19  m/sec. 

Other  Complications.  In  addition  to  the  gust-scale,  instrumental, 
advective,  and  non-steady  state  wind  variations,  two  other  complications 
arise  in  the  planetary  boundary  layer  to  prevent  the  observation  of  the 
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Fig.  9.  Normalized  frequency  distributions  of  lapse  rate  in  surface 
900  mb  layer  (from  Mendenhall,  1967). 

simple  Ekman  profile.  These  are: 

1.  Vertical  stability  variations 

2.  Horizontal  temperature  gradient  variations  '^dermal  wind  or 
geostrophic  veering  influence) 

Stability  Variations.  Stability  influences  are  a  major  factor  in  the 
veering  character  over  land.  Over  the  oceans,  diurnal  stability  varia- 
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Fig.  10.  Frictional  veering  (with  geostrophic  veering  effect  accounted 
for  and  removed)  as  a  function  of  lapse  rate  in  the  surface  to 
900  mb  layer  at  each  station.  Smooth  lines  were  obtained  by 
averaging  frictional  veering  in  classes  of  lapse  rate  and  con¬ 
necting  the  points  (from  Mendenhall,  1967). 

tions  are  small  and  do  not  significantly  influence  the  veering.  Fig.  9 
shows  the  much  larger  lapse  rats  variations  at  Shreveport  (a  typical 
land  station)  compared  with  xapse  rate  variations  at  weather  ships  I,  J, 
N  and  Johnston  Island.  Fig.  10  shows  the  variations  of  frictional  veer¬ 
ing  with  lapse  rate  in  the  surface -to- 900  mb  layer  for  these  same  sta¬ 
tions.  Veering  is  directly  related  to  stability.  On  the  average  the  1st 
km  oceanic  veering  decreases  about  10°  for  stability  changes  from  iso¬ 
thermal  to  dry-adiabatic. 

Stability  variations  over  the  oceans  are  typically  small.  The  influ¬ 
ence  of  p.  b.  1.  oceanic  lapse  rate  variations  are  thus  not  felt  to  be  a 
dominant  influence  in  determing  differences  in  observed  veering. 
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Fig.  11.  Observed  veering  in  the  surface  to  900  mb  layer  at  Ship  N  as 
a  function  of  geostrophic  veering,  obtained  by  averaging  ob¬ 
served  veering  in  5°  class  intervals  of  geostrophic  veering. 
The  heavy  line  connects  these  average  values.  The  dashed 
lines  indicate  that  the  median  ("50%")  geostrophic  veering  is 
-8°  and  that  the  median  observed  veering  is  2°.  The  differ¬ 
ence  of  10°  is  the  median  frictional  veering.  The  thin  sloping 
line  connects  points  of  equal  observed  and  geostrophic  veer¬ 
ing,  or  where  frictional  veering  equals  zero.  The  vertical 
distance  between  the  lines  thus  represents  frictional  veering 
and  is  indicated  by  arrows  at  the  median  point.  The  "75%" 
mark  indicates  that  75%  of  the  observations  had  geostrophic 
veering  of  less  than  3°,  etc.  Sample  size:  2386  observations 
(from  Mendenhall,  1967). 


Horizontal  Temperature  Gradient  Variations.  In  contrast  the  ther¬ 
mal  wind  influence  (or  geostrophic  veering)  can  be  very  pronounced  with 
cold  and  warm  air  advection  and  must  always  be  accounted  for.*  Ther¬ 
mal  wind  or  geostrophic  veering  influences  are  felt  to  be  largely  elim- 


It  is  only  the  temperature  gradient  along  the  direction  of  motion 
which  significantly  alters  the  measured  veering. 


Q 

£ 

CE 

8 

®  -10 


-10  -6  0  10  20 
GEOSTROPHIC  VEERING  (deg.) 

Fig.  12.  Observed  veering  in  the  surface  to  900  mb  layer  at  Johnston 
Island  as  a  function  of  geostrophic  veering.  (See  Fig.  11  for 
further  description).  Dashed  lines  indicate  median  geostro¬ 
phic  veering  of  -6°  and  median  observed  veering  of  4°.  Dif¬ 
ference  of  10°  is  median  frictional  veering.  Sample  size: 
3667  observations  (from  Mendenhall,  1967). 


Fig.  13. 


*'0  o  10  20  30 

GEOSTROPHIC  VEERING  (deg.) 

Observed  veering  in  the  surface  to  900  mb  layer  at  Ships  I  and 
J  asa  function  of  geostrophic  veering.  (See  Fig.  11  for  furth¬ 
er  description).  The  median  frictional  veering  of  7°  is  in¬ 
dicated  bv  arrows.  (Original  data  after  Findlater  et  al. ,  1966; 
from  Mendenhall,  1967). 
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inated  by  assuming  constant  baroclinicity  '  through  Uic  lowest  two  kilo¬ 
meters  and  then  subtracting  die  veering  of  the  second  km  being  assumed 
to  be  a  product  of  the  thermal  wind  alone.  If  the*  observed  wind  veering 
in  the  first  and  second  km  layer  is  I  f>  and  6°,  respectively,  then  the 
thermal  wind  influence  would  be  assumed  to  lie  6°.  This  a°  veering 
would  be  subtracted  from  the  first  km  veering  data  and  a  frictional  veer¬ 
ing  of  10°  would  lie  accepted.  This  procedure  has  been  verified  by  mea¬ 
suring  the  1st  km  horizontal  temperature  gradients  directly  and  deter¬ 
mining  how  this  gradient  compares  with  the  one  obtained  by  subtracting 
the  second  km  veering.  It  has  been  found  to  be  generally  satisfactory 
but  not  representative  at  stations  with  an  inversion  layer  below  800  mb. 

The  subtraction  of  the  thermal  wind  influence  can  reduce  and  norm¬ 
alize  the  frictional  veering  as  has  been  demonstrated  by  Mendenhall, 

1967  (see  Figs.  1  1-18).  Fig.  11  shows  the  frequency  distribution  of  ob¬ 
served  veering  as  a  function  of  geostrophic  veering  (or  thermal  wind) 
at  weather  ship  N.  Similar  plots  of  observed  vs.  geostrophic  veering 
are  shown  in  Figs.  12  and  18  for  Johnston  Island  and  weather  ships  I 
and  J  (data  from  Findlater  et  al.  ,  1966).  Note  how  large  the  thermal 
wind  correction  can  be.  Note  also,  that  after  it  is  made,  the  average 
frictional  veering  is  10°. 

The  small  net  veering  at  Ship  N  in  Fig.  1  is  due  to  cold  air  advec- 
tion.  Warm  air  advection  is  occuring  at  Swan  Island  (Fig.  2).  Subtrac- 

.j.* 

Thermal  wind  influence,  geostrophic  veering,  baroclinicity  are 
synonymously  used. 
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tion  of  these  opposite  acting  thermal  wind  influences  largely  reduces 
these  veering  differences. 

Summary.  Observed  wind  veering  is  a  function  of  a  number  of  sim¬ 
ultaneously  acting  influences.  Thus, 


Observed  Frictional  Geostrophic  + 
Veering  Veering  or  Thermal 

Wind  Veer¬ 
ing 

A  =  B  +  C  + 


Stability+  Unrepre-+  Local  and 
Change  sentative  Advective 
Veering  Turbulent  Veering 
and  Instru-  Change 
ment  Veer¬ 
ing 

D  +  E  +  F 


The  attempt  of  this  paper  is  to  determine  B  i  om  large  statistical 
samples  of  A.  E  veering  influences  are  eliminated  by  averaging  over 
a  large  data  sample.  D  and  F  veering  influences  are  accounted  for  or 
eliminated  by  averaging  over  a  large  data  sample,  a  long  period,  and  a 
large  longitude  (360°)  belt.  C  veering  influences  are  eliminated  or  ac¬ 
counted  for  by  large  time  and  space  averaging  and  by  assuming  constant 
baroclinicity  with  height  and  subtracting  the  2nd  km  veering  from  the  1st 
km.  This  is  shown  to  be  very  effective.  Thus,  B  is  determined  from 
A  after  the  veering  influences  of  C,  D,  E,  and  F  have  been  eliminated 


or  accounted  for  and  subtracted  out. 


III.  DATA  SAMPLE 


Two  primary  oceanic  upper  air  data  sources  are  involved:  islands 
and  atolls,  and  ships. 

Islands  and  Atolls  Data.  This  data  consists  of  regular  island  or  atoll 
upper  air  reporting  sites  from  which  two  or  four  daily  soundings  are  a- 
vailable.  This  stationary  group  of  observations  has  been  analyzed  in 
three  groups. 

1)  fifteen  of  the  stations  were  sampled  continuously  for  a  6  month 
summer  period  in  the  early  1960's.  Table  1  summarizes  the 
imormation  on  these  15  stations  obtained  from  the  U.S.  Weather 
Bureau  Northern  Hemisphere  Data  Tabulations. 

2)  veering  information  was  extracted  from  tropical  stations  during 
1967-69  in  the  Western  Pacific  and  in  the  West  Indies-Gulf  of 
Mexico  region.  This  data  v  as  obtained  in  connection  with  the 
large  statistical  studies  of  tropical  cloud  clusters  (Williams, 

1970)  and  clear  areas  which  the  author  and  his  graduate  students 
have  been  investigating.  The  data  was  stratified  only  by  ocean 
and  latitude  (greater  than  or  less  than  18°).  Wind  information 
was  obtained  from  magnetic  tape  images  of  card  deck  645  of  the 
National  Climatic  Center,  Asheville,  North  Carolina,  at  the 
surface,  1000,  95 J,  900,  850  and  800  mb.  Winds  are  reported 
to  the  nearest  m/sec  and  in  direction  to  the  nearest  degree. 

3)  a  smaller  portion  of  the  1967-69  data  was  separated  into  indivi¬ 
dual  cloud  cluster  and  clear  area  veering  groups.  Data  source 
is  the  same  as  in  2).  Figs.  14  and  15  show  the  locations  of  the 
stations  used  in  the  latter  two  data  groups. 

80,000  Ocean  Ship  Veering  Reports  from  1949-1964.  Nearly  all 
U.S.  National  Climatic  Center,  Asheville,  North  Carolina  data  on  the 
variation  of  wind  in  the  lowest  two  km  over  the  oceans  has  been  assem¬ 
bled  for  the  16-year  period  of  1949-64.  Data  includes  over  80,000  ship 
rawin  and  pibal  observations  from  merchant  ships  which  took  rawins  or 
pibals,  weather  ships,  military  ships--every  kind  of  surface  vessel 
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Table  1 

Geographical  Data  on  Stationary  Locations 


Name  Location  Months  of  Highest  Height  (meters)  and 

Observations  Point  Location  of  Observing 
(meters)  Station 


I'ac 

i  f  i  c 

Marcus  I . 

24N , 

l  Vi  K 

6 

22 

20, 

center  part 

Midway  I . 

28N, 

1  77W 

f> 

14 

13, 

center 

French  Frigate  Shoals 

24N, 

1 6(>W 

6 

(volcanic 

rock) 

3 

T-hnston  I. 

16N, 

1 70W 

6 

22 

2, 

north  side 

Wake  I. 

20N , 

166E 

5 

7 

4, 

SE  side 

Eniwetok  Atoll 

1  IN, 

162E 

6 

4 

3, 

SE  side 

Kwajalein  Atoll 

9N, 

168E 

6 

3 

3, 

center 

Majuro  Atoll 

;n. 

171E 

6 

5 

5, 

west  side 

Palmyra 

5N, 

162W 

1% 

3 

2 

Weather  Ship  V 

34N, 

164E 

6 

25 

10 

Atlantic 


Grand  Turk  I. 

22N,  71W 

6 

10 

5 

Grand  Cayman  I. 

19N,  81W 

6 

10 

5 

Swan  1. 

17N,  84W 

6 

10 

5 

San  Andrds  I. 

12N,  82W 

6 

110 

5 

Weather  Ship  E 

35N,  48W 

10 

25 

10 

21 


Fig.  15.  Stations  from  which  veering  information- in  Tables 

5  and  6  was  obtained. 


which  Asheville  had  on  record. 

Between  1949  and  1955  (data  set  I)  wind  reports  were  available  for 
the  surface,  500  m,  1  km,  1.  5  km  and  2  km  and  directions  were  report¬ 
ed  on  the  16  point  compass  (22-1/2°  intervals).  A  little  less  than  half 
the  data  is  in  this  category.  The  reports  of  the  1950  through  1964  per¬ 
iod  (data  set  II)  give  wind  information  to  the  nearest  degree  and  by 
height  at  the  surface,  150,  300,  500,  1  km,  1.5  km  to  2  km.  Veering 
was  obtained  from  the  earlier  (1949-55)  16  point  compass  data  by  as¬ 
suming  the  mean  wind  direction  of  each  16  point  direction  category  and 
noting  changes  of  category  with  height.  In  the  statistical  average  this 
proved  very  satisfactory. 

Data  stratifications  were  performed  on  each  of  the  two  data  sets  sep 


arately.  At  all  levels  and  in  all  stratifications,  only  very  small  data 


sot  differences  were  found.  To  increase  the  sample  size,  the  author 
felt  quite  justified  in  averaging  both  data  samples  together.  This  was 
accomplished  l  y  500  m  height  intervals.  Figs.  16-18  give  the  number 
of  observations  available  in  each  1 0°latitude- longitude  area.  These  fig¬ 
ures  show  how  tin  two  ship  data  sets  and  the  combined  data  sample  (sets 
I  and  II)  are  distributed  geographically.  Data  is  concentrated  along  the 
shipping  routes,  but  a  fairly  good  latitudinal  and  longitudinal  distribu¬ 
tion  is  present. 

Some  of  the  ship  reports  had  to  be  eliminated  because  of  missing  in¬ 
formation  at  some  levels,  or  if  speeds,  shears,  or  veering  angles  were 

unrealistically  excessive.  After  this  screening  analysis  there  were 
60,000  ship  reports  remaining.  Over  5000  observations  are  available 

in  the  latitude  Ire  It  from  10°-20°N.  This  ship  information  has  been 
stratified  by  latitude,  season,  wind  direction  and  wind  speed.  Compon¬ 
ent  wind  shears  along  and  perpendicular  to  the  wind  at  2  km  have  been 
printed  out  by  :>00  m  intervals.  In  addition,  averaging  was  performed 
to  obtain  the  wind  speed  squared  and  cubed.  Vertical  gradients  of  this 
information  were  taken  at  the  various  reporting  levels.  Fig.  19  shows 
a  typical  data  printout  of  this  ship  vessel  information  for  data  set  II. 
Symbols  are  defined  in  Table  2  and  the  stratifications  which  were  made 
are  shown  in  Table  3. 

Since  1964  this  class  of  ship  data  is  available  only  at  50  mb  inter¬ 
vals.  This  type  of  vertical  resolution  is  no  longer  available. 
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(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ID 


Table  2 

DEFINITION  OF  SYMBOLS 

£Vi 

SB  -  The  average  value  of  the  wind  speed;  — jj-  ;  [Units-m/sec]; 
N  «  #  of  obs.;  V  =  speed 

52  -  The  average  value  of  the  square  of  the  wind  speed; 

[Units-m1 2 3/sec2] 

£(Vi)3 

53  -  The  average  value  of  the  cube  of  the  wind  speed;  — ^ 

[Units-m^/sec2] 


-  The  component  of  the  resultant  vector  wind  along  the  x-axis 

I  i 

(in  the  meteorological  coordinate  system).  — jj-  [Units-m/sec] 

-  The  component  of  the  resultant  vector  wind  along  the  y-axis 

£vi 

(in  the  meteorological  coordinate  system).  — jj- 


UA  -  The  average  shear  (between  the  indicated  levels)  parallel  to 
the  wind  direction  at  the  lower  level.  A  positive  UA 
indicates  the  component  of  the  wind  along  the  direction 
of  flow  at  the  lower  level  is  increasing  with  height. 


VA  -  The  average  shear  (between  the  indicated  levels)  perpendicular 
to  the  wind  vector  at  the  lower  level.  Here  a  positive 
VA  indicates  the  wind  is  veering  with  height. 


VS  -  The  average  of  the  absolute  values  of  the  observed  vertical 
wind  shear ;  „  .  u  -  (VA)2 

ii 

DD  -  The  wind  veering  computed  from  the  values  of  VA  and  SB;  generally 
VA 

DD  »  ATAN(— )  x  57  3  Here  the  value  of  SB  at  the 

lower  level  was  used. 

DA  -  Mean  wind  veering  computed  using  wind  direction  data  and  the 
following  limits: 


(1)  All  calm  winds  were  neglected. 

(2)  All  cases  with  winds  greater  than  30  m/sec  at  any 
level  were  neglected. 

(3)  All  cases  with  direction  changes  greater  than  45° 
between  two  successive  levels  or  greater  than  90° 
between  any  of  the  levels  were  neglected. 


DR  -  Mean  wind  veering  computed  using  wind  direction  data  and  the 
following  limits: 

(1)  All  winds  less  than  3  .n/sec  were  neglected. 

(2)  Other  criteria  are  the  same  as  those  in  (2)  and  (3)  for  DA. 


period  1956-1964.  Symbols  defined  in  Table  2. 
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Preceding  page  blank 

Table  3 

. 

STRATIFICATIONS 

5 

1. 

ALL  WINDS  BY  LATITUDE 

! 

2. 

NORTH  WINDS  BY  LATITUDE 

3. 

EAST  WINDS  BY  LATITUDE 

4. 

SOUTH  WINDS  BY  LATITUDE 

1 

5. 

WEST  WINDS  BY  LATITUDE 

6. 

WINTER  WINDS  BY  LATITUDE 

7. 

SPRING  WINDS  BY  LATITUDE 

8. 

SUMMER  WINDS  BY  LATITUDE 

i 

9. 

FALL  WINDS  BY  LATITUDE 

- 

10. 

<  9  m/s ec  WINDS  BY  LATITUDE 

11. 

~  9  m/sec  WINDS  BY  LATITUDE 

' 

12. 

0-4  ra/sec  WINDS  BY  LATITUDE 

13. 

5-7  m/sec  WINDS  BY  LATITUDE 

14. 

8-11  m/sec  WINDS  BY  LATITUDE 

15. 

~  12  m/sec  WINDS  BY  LATITUDE 

; 

Wind  direction  classification  (based  on  direction  at  500m  level): 

Set  I.  N  -  (16,  1,  2)  Set  II.  N  -  316°-45° 

__  E  -  (4,  5,  6)  1Q,,  E  -  46°-135°  .  _  . 

1949-55  v  1956-o4  points  of  compass 

S  -  (8,  9,  10)  S  -  136°-225° 

W  -  (12,  13,  14)  W  -  226°-315°  classification 

Season  classification: 

• 

Winter  --  January,  February,  March 

Spring  —  April,  May,  June 

Summer  —  July,  August,  September 

Fall  —  October,  November,  December 

. 

A 
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IV.  WIND  VEERING  AT  ISLANDS  AND  ATOLLS 

Individual  Station  Long  Term  Veering.  Table  4  lists  the  long  term 
(mostly  summer  months  except  for  the  one  and  half  month  average  at 
Palmyra  during  the  Line  Island  Experiment)  wind  veering  in  the  1st  and 
2nd  km  layers  for  15  individual  island  and  atoll  stations  in  the  tropics 
and  subtropics  as  described  in  Table  1.  Observations  were  taken  every 
6  or  12  hours.  More  than  350  observations  are  available  for  most  sta¬ 
tions.  Figs.  20  and  21  show  the  6  morth  average  vertical  profiles  of 
the  variation  of  angle  veering  and  wind  speed  through  the  lowest  2  km 
layers.  Note  the  substantial  differences  between  the  1st  km  veering  (12° 
average)  and  the  2nd  km  veering  (1°  average).  If  the  frictional  veering 
is  assumed  to  be  approximated  by  the  difference  between  the  1st  and  2nd 
km  layer  veering,  then  a  frictional  veering  of  about  11°  results.  The 
assumption  of  constant  thermal  wind  and  the  subtracting  of  the  2nd  km 
veering  from  the  1st  helps  to  reduce  the  1  km  veering  spread.  This  is 
particularly  so  with  the  individual  monthly  veering  values. 

Individual  Station  Mean  Monthly  Veering.  Figs.  22-29  portray  the 
individual  monthly  average  angle  veerings  for  8  of  these  15  stations. 

Note  how  the  individual  monthly  veerings  can  be  significantly  different 
from  their  6  month  averages.  Difference  in  monthly  thermal  wind  (as 
sea-surface  temperature  gradients  change)  are  felt  to  be  primarily  re¬ 
sponsible  for  these  monthly  veering  differences.  Lapse  rate  variations 
may  also  occur,  but  are  believed  to  be  only  a  secondary  cause  for  these 
differences.  It  can  be  seen  that  subtraction  of  the  monthly  average  2nd 
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Table  4 

Veering  Data  at  Stationary  Locations 
Average  Veering  (degrees) 


Name 


Months  of  (1  km-sfc)  (2  km-1  km)  (1  km-sfc)- 

Observations  _  _  (2  km-1  km) 

(2  or  4  obs.  per  day) 


Marcus  I. 

6 

Pacific 

11 

0 

11 

Midway  I. 

6 

16 

-2 

18 

French  Frigate  Shoals 

6 

3 

-4 

7 

Johnston  I. 

6 

5 

-2 

7 

Wake  I. 

5 

13 

0 

13 

Eniwetok  Atoll 

6 

4 

-1 

5 

Kwajalein  Atoll 

6 

12 

1 

11 

Majuro  Atoll 

6 

12 

2 

10 

Weather  Ship  V 

6 

15 

2 

13 

Palmyra 

1*5 

9 

0 

9 

Grand  Turk  I. 

6 

Atlantic 

17 

2 

15 

Grand  Cayman  1. 

6 

15 

3 

12 

Swan  I. 

6 

21 

2 

19 

San  Andr€s  1. 

6 

12 

4 

8 

Weather  Ship  E 


10 


9 


1 


8 


WIND  DIRECT  DIFFERENCE  RELATIVE  TO  2 Km  HEIGHT 


T  ig.  20.  Six  month  average  of  individual  station  wind  veering 
Solid  line  represents  the  average  of  the  15  stations. 


HEIGHT 
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VELOCITY  (m/sec) 


Fig.  21.  Six  month  average  of  individual  station  tangential  wind 
Solid  line  represents  the  average  of  the  15  stations. 
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km  veering  from  the  1st  km  veering  does  help  to  reduce  the  spread  of 

the  profiles. 

Veering  Information  from  West  Pacific  and  West  Indies.  This  data 
was  obtained  in  connection  with  the  author's  tropical  cloud  cluster  and 
clear  area  studies  from  stations  shown  in  Figs.  14  and  15.  The  data  is 
not  continuous  and  the  stations  have  been  mixed.  Roughly  equal  amounts 
of  data  are  available  from  all  stations.  Table  5  lists  the  1st  and  2nd  km 
mean  veering  information  for  these  two  regions.  The  West  Indies  data 
has  been  divided  into  two  latitude  belts  (greater  and  less  than  18°).  It 
is  again  to  be  noted  that  the  1st  and  2nd  km  veerings  average  about  12° 
and  2°.  The  West  Indies  data  north  of  18°  shows  larger  veerings. 

Fig.  30  is  the  typical  vertical  profile  of  the  veering  from  West  Pa¬ 
cific  atolls  and  selective  surface  ships  surrounding  the  atolls.  This  fig¬ 
ure  shows  that  there  is  essentially  no  difference  between  the  statistical 
averages  of  the  island-atoll  and  the  ship  veering. 

Cloud  Cluster  vs.  Clear  Area  Veering.  Table  6  shows  the  1st  and 
2nd  km  wind  veering  for  stations  within  2°  latitude  of  the  center  of  sat¬ 
ellite-observed  tropical  cloud  clusters,  and  the  same  veering  information 
relative  to  tropical  clear  regions.  Note  that  for  the  clusters  the  1st  km 
veering  is  more  (15°  vs.  11°),  but  it  extends  well  into  the  second  km 
level,  (6°  veering  for  the  clusters  vs.  2°  veering  for  the  clear  areas). 
This  is  to  be  expected  if  the  cumulus  clouds  act  to  carry  the  turbulent 
sub-cloud  layer  momentum  to  higher  levels.  The  clear  area  veering  is 
very  similar  to  the  average  of  the  ship  and  atoll-island  stations.  When 


Table  5 


1st  and  2nd  km  VEERING  FROM  TROPICAL  AND  SUBTROPICAL 


ATOLL  AND  ISLAND 

STATIONS 

1  km-sfc 

2  km-1  km 

(1  km-sfc)  - 
(2  km-1  km) 

Western  Pacific 

7210  cases 

10° 

2° 

8° 

West  Indies  (less  18° 
9205  cases 

latitude)  10° 

3° 

7° 

West  Indies  (greater 
7980  cases 

18°  latitude)  14° 

3° 

11° 

AVERAGE 

11° 

3° 

8° 

Table 

6 

CLOUD  CLUSTER  vs.  CLEAR 

AREA 

1st 

and  2nd  kn 

i  VEERING 

CLUSTERS 

(1  km- 

sf  c) 

(2  km-1 

.  \  (1  km-sfc) 
k")  (2  km-1  km) 

Western  Pacific 
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Fig.  30. 

Observed  wind  angle  veering 
with  height  in  the  lowest  two 
km  for  atoll  and  ship  reports 
of  the  tropical  Pacific.  Wind 
direction  at  two  km  is  used  as 
reference.  Curve  a  repre¬ 
sents  veering  of  wind  with  height 
from  ships  which  were  located 
at  least  1°  latitude  from  any 
land.  Curve  b  represents  the 
frictional  veering  of  wind  with 
height  as  observed  from  atoll 
data;  curve  c  as  observed 
from  ships  located  within  1° 
latitude  of  land.  Number  of 
observations  in  each  class  is 
shown  on  the  figures. 


OC  (Degrees) 


the  2nd  km  veering  is  subtracted  from  the  1st  km  veering  an  equal  9° 
veering  is  observed  for  both  the  clusters  and  the  clear  areas. 

Conclusion  -  The  long  term  average  from  the  tropical  and  sub-trop¬ 


ical  islands  and  atolls  indicate  a  typical  1st  km  'frictional  veering'  as 
here  defined  of  about  10°. 
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V.  ANGLE  VEERING  FROM  OCEAN  SHIPS 

Veering  vs.  Latitude.  Figs.  31  and  32  portray  the  latitude  distri¬ 
bution  of  the  observed  wind  angle  veerings  for  the  lowest  two  km  layers 
for  data  sets  I  (1949-55)  and  II  (1956-64).  There  appears  to  be  no  signi¬ 
ficant  latitude  difference  in  the  veerings  of  either  data  sample.  The  1st 
km  layer  shows  a  distinct  10°  greater  veering  than  the  2nd  km  layer. 
Figs.  33  and  34  show  the  latitude  distribution  of  the  surface  (10  m  ship 
level)  and  1  km  wind  components  perpendicular  to  the  two  km  wind. 
Positive  values  indicate  veering.  The  one  and  two  km  winds  have  very 
small  perpendicular  component  difference.  The  surface  component  per¬ 
pendicular  to  the  2  km  wind  averages  about  1  m/sec.  The  p.b.l.  "fric¬ 
tional  veering11  due  to  mechanically  disturbed  gust-scale  motions  is  now 
defined  as  the  difference  between  the  1st  and  2nd  km  layer  veering  an¬ 
gles  as  given  in  Figs.  31-32.  This  veering  is  also  proportional  to  the 
perpendicular  components  as  seen  by  the  hatched  area  in  Figs.  33  and 
34. 

Observing  that  there  is  very  little  difference  in  the  data  samples  (as 
previously  mentioned)  the  data  of  both  samples  were  combined.  Fig.  35 
shows  angle  veering  by  500  m  layers  from  the  2  km  reference  level 
down  to  the  surface.  The  small  difference  in  veering  angle  variations 
with  latitude  and  the  much  larger  1st  km  layer  veering  are  clearly  seen. 
Fig.  36  graphically  portrays  the  entire  oceanic  Hemisphere  average  of 
the  1st  and  2nd  km  layer  veering  angle  differences.  These  average 


Fig.  31.  Latitude  variation  of  the  average  rawin  and  pibal  ship 
wind  veering  in  first  and  second  km  layer  for  1949- 
1955  data  where  wind  direction  was  reported  to  the 
16  point  compass. 


LATITUDE 


Fig.  32.  Latitude  variation  of  the  average  rawin  and  pibal  ship 
wind  veering  in  the  first  and  second  km  layer  for  data 
of  1956-1964  where  wind  direction  was  reported  to  the 
nearest  degree. 
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LATITUDE 


Fig.  33.  Latitude  distribution  of  surface  (sfc)  and  1  km  wind 
components  (  m/sec  )  perpendicular  to  the  second 
km  wind  for  ship  vessel  data  of  1949-1955.  Hatched 
area  portrays  the  assumed  frictional  veering  which 
is  the  difference  in  surface  and  1  km  components. 


LATITUDE 


Fig.  34.  Latitude  distribution  of  surface  (sfc)  and  1  km  wind  com 
ponent  (m/sec)  perpendicular  to  the  second  km  wind  for 
ship  vessel  data  of  1956-1P64.  Hatched  area  portrays 
the  assumed  frictional  veering  which  is  the  difference 
in  the  surface  and  1  km  components. 
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VEERING  DEGREES  FROM  2km  WIND 


NO.  OF  OBS. 

J. - |  5166  |  8011  |  24.390  j  17.436  j  4432  |  324 

Eq.  10°  20°  30°  40°  50°  60° 

LATITUDE 


Fig.  35.  Latitude  distribution  of  wind  direction  veering  in  500  m 
layers  relative  to  the  2  km  level  wind  for  the  entire  data 
sample  of  1949-1964. 


WIND  VEERING  IN  FIRST 
AND  SECOND  KM  LAYER 


Fig.  36.  Average  of  first  and  second  km  wind 
veering  for  all  ship  vessel  data. 


O*  11-1/2°  and  2-1/2°  giving  a  mean  oceanic  frictional  veering  in 
the  1st  km  layer  of  ^  9°. 

Veering  vs.  Wind  Direction  and  Latitude.  Fig.  37  portrays  the  large 
difference  of  1st  km  wind  veering  when  the  reports  are  stratified  by  wind 
direction  to  the  four  point  compass  as  shown  in  Fig.  38.  Note  the  15° 
to  25°  1st  km  veering  differences  between  north  and  south  wind,  yet 
east  and  west  winds  show  only  small  veering  differences.  These  north 
and  south  wind  veering  differences,  believed  to  be  due  to  thermal  wind 
difference  in  cases  of  warm  and  cold  air  advection,  are  mostly  elim¬ 
inated  when  the  2nd  km  veering  is  subtracted  from  the  1st  km  wind 
with  constant  baroclinicity  assumption.  Fig.  39  shows  both  the  1st  and 
2nd  km  wind  veering  for  north  and  south  winds.  When  the  2nd  km  veer¬ 
ing  is  subtracted  from  the  1st  km  veering  the  now  defined  'frictional 
veering'  is  the  same  for  both  directions.  The  20°  and  10°  average  1st 
and  2nd  km  veerings  with  south  winds  give  a  similar  frictional  veering 
as  the  0°  and  -10°  1st  and  2nd  km  veerings  with  north  winds.  Table  7 
portrays  this  veering  information  in  20°  latitude  segments. 

It  is  thus  seen  that  there  is  no  observed  major  1st  km  'frictional 
veering'  difference  with  wind  direction  or  latitude. 

Veering  vs.  Wind  Speed  and  Latitude.  Figs.  40-43  show  the  1st  and 
2nd  km  veering  difference  (ie.  frictional  veering)  by  latitude  for  wind 
speed  categories  of  0-4,  5-7,  8-11  and  greater  than  12  m/sec.  It  is 
observed  that  there  are  no  primary  'frictional  veering'  differences  as¬ 
sociated  with  the  four  speed  categories  at  ail  latitudes.  Fig.  44  shows 
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COMBINED  RAWIN  AND  PJBAL  OBS. 

OBSERVED  WIND  VEERING  (DEGREES)  IN  Sfc  *IKm 
LAYER  BY  WIND  DIRECTION  AND  LATITUDE 


EQ.  10*  20*  30*  40»  50*  60»  P 

LATITUDE 

Fig.  37. 


Fig.  38.  Wind  direction  stratification. 
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NORTH  AND  SOUTH  WIND  VEERINGS  (Degrees) 
1ST  AND  2ND  Km  LAYERS  AND  DIFFERENCE  - 
assumed  to  be  the  frictional  veering. 


i—  32* 


|— 24° 


-  16° 

-  8* 

-  O* 

--8* 

—16* 


EQ.  10“  20®  30°  40“  50°  60“ 

LATITUDE 


Fig.  39.  Portrayal  of  1st  and  2nd  km  wind  veering 
differences  for  north  and  south  winds. 

that  there  is  no  frictional  veering  difference  with  wind  speed  in  the 
tropics  and  subtropics.  In  the  westerly  latitudes,  however,  strong 
winds  show  a  frictional  veering  about  twice  that  of  the  weak  winds.  The 
magnitude  of  the  2nd  km  veering  is  larger  for  higher  wind  speeds  as 
was  expected  (See  Discussion  Section).  Table  8  more  explicitly  shows 
this  veering  stratification  by  wind  speed. 


WIND  VEERING  (DEGREES)  „  WIND  VEERING  (DEGREES) 


43 


20*  r- 


16°  (— 


Combined  Rawin  8  Pibal  Obs. 

Wind  Veering  in  First  and  Second  km.  Layer  - 
-For  Wind  Speeds  0-4  m/sec. 


12® 


_  Sfe  to  I  km 


8®  h 


4®  h- 


-4® 


112  i  672  i 


NO.  OF  OBS 
1540  i  4911  i  1579  i  286  I 


EQ 
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30® 

LATITUDE 


40°  50“  60“ 


^g.  40.  Latitude  distribution  of  frictional  veering  (hatched 
area)  for  500  m  winds  of  0-4  m/sec. 


Fig.  41.  Latitude  distribution  of  frictional  veering  (hatched 
area)  for  500  m  winds  of  5-7  m/sec. 
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LATITUDE 


Fig.  42.  Latitude  distribution  of  frictional  veering  (hatched 
area)  for  500  m  winds  of  8-11  m/sec. 


LATITUDE 


Fig.  43.  Latitude  distribution  of  frictional  veering  for  300  m 
winds  >  12  ni/sec. 
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Preceding  page  blank 


WIND  VEERING  (Deg.)  by 
WIND  DIRECTION  and  LATITUDE 


Eq-  20° 

20°  -  40* 

40°-  60* 

global 

MEAN 

from  North 

8.3 

1  km  -  Sfc. 

-1.1 

1.9 

-  3 

South 

14.1 

20.3 

26.5 

-  20 

East 

9.4 

11.2 

17.5 

-  13 

West 

8.0 

7.1 

7.9 

-  8 

Annual 

10.3 

k7 

12.6 

-  lO'j 

from  North 

-0.8 

2  km  -  1  km 

-10.6 

-7.4 

~  -6 

South 

5.6 

12.1 

1 1.0 

-  10 

East 

3.8 

3.0 

9.8 

-  6 

West 

-0.3 

l.l 

0.3 

-  0 

Annual 

2.9 

1.2 

2.3 

-  H 

Frictional  Veerin 

g  ( 1  Km  -  S 

fc.)  -(2km 

-  Ikm) 

from  North 

7.5 

11.7 

93 

-  9 

South 

8.5 

8.2 

15.5 

-10 

Eost 

5.6 

8.2 

7.7 

-  7 

West 

8.3 

6.0 

7.6 

^  8 

Annual 

7.4 

7.5 

10.3 

-  8 

Table  7. 

Veering  vs.  Season  and  Latitude.  Figs.  45-48  show  the  latitude 
variation  of  1st  and  2nd  km  veering  by  season  and  Fig.  49  portrays  1st 
km  'frictional  veering'  in  combination  for  all  seasons.  Table  9  more 
explicitly  lists  these  seasonal  veering  values  by  20°  latitude  intervals. 
Table  10  shows  the  1st  and  2nd  km  perpendicular  component  variations 
of  veering  by  season. 
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Fig.  44.  Latitude  distribution  of  frictional  veering  by  wind  speed. 

Seasonal  differences  are,  in  general,  not  very  large.  The  lower 
autumn  and  winter  frictional  veering  in  middle  latitudes  are  felt  to  be 
partly  related  to  increased  vertical  lapse  rate  instability  brought  about 
by  cold  air  moving  over  warmer  water  in  the  shipping  lanes  of  the  west¬ 
ern  oceans  where  the  ship  data  is  concentrated.  Fig.  10  shows  that 
one  would  expect  a  decrease  of  observed  veering  with  increasing  vert¬ 
ical  instability. 

It  is  thus  seen  that  there  are  no  large  seasonal  veering  differences. 

Conclusion.  Ship  observations  indicate  that  the  1st  km  oceanic  'fric¬ 
tional  veering'  has  no  primary  relationship  with  latitude,  wind  direction. 
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WIND  VEERING  (Deg.)  by 
WIND  SPEED  and  LATITUDE 


WIND  SPEEDS 

m/sec 

Eq  -  20 
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-  2 

3-11 

3.2 

-0.2 

1.8 

-  3 

212 

4.4 

2.4 

1.0 

-  4 

Annual 

2.9 

1.2 

2.3 

-  3 
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b  6 
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1 1.0 

-  9 
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5.8 

82 
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7.5 

10.3 
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Table  8. 


season,  or  with  wind  speed  in  tropical  and  subtropical  latitudes.  In 
westerly  latitudes,  strong  winds  veer  more  than  weak  winds.  For  all 
stratifications,  the  first  km  observed  veering  is  observed  to  be  about 
10-11°,  the  second  km  veering  about  2°.  These  results  hold  very  well 
even  at  deep  tropical  latitudes.  The  lack  of  any  observed  frictional 
wind  veering  correlation  with  latitude  is  most  striking  and  significant. 
This  is  discussed  in  the  last  section. 


WIND  VEERING  (DEGREES; 
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Combined  Rawin  8  Pibol  Obi 
Observed  Wind  Veering  in  First  and 
Second  km.  Lover  -  Winter 


Sfc.  to  I  km. 


Frictional 

Veering 


A  to  2  km 


NO  OF  OBS 


20°  30° 

LATITUDE 


40“  50°  60° 


Fig.  45.  Latitude  distribution  of  frictional  veering  in  winter. 


Combined  Rawin  8  Pibal  Obs. 
Observed  Wind  Veering  in  First  and 
Second  km  Layer  -  Spring  y 
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LATITUDE 


Fig.  46.  Latitude  distribution  of  frictional  veering  in  spring. 


53 


LATITUDE 


Fig.  47.  Latitude  distribution  of  frictional  veering  in  summer, 


LATITUDE 

Fig.  48.  Latitude  distribution  of  frictional  veering  in  autumn. 


WIND  VEERING  (Deg.)  by 
SEASON  and  LATITUDE 
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Table  9. 


Combined  Rawin  8  Pibal  Obs. 
Frictional  Veering  (Degrees)  by 
Season  and  Latitude 
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Flff  49.  Frictional  veering  by  latitude  and  season, 
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VI.  DYNAMIC  CONSIDERATIONS 

Wind  Speeds.  The  average  speed  of  the  ship  surface  winds  (uQ)  as 
measured  at  10  meters  (deck  level)  height  are  shown  in  Fig.  50.  Typ¬ 
ical  values  are  6-8  m/sec,  less  in  deep  tropics,  more  in  westerly  lat¬ 
itudes.  The  surface  stress  and  kinetic  energy  dissipation  are  roughly 
proportional  to  the  square  and  cube  of  the  surface  winds,  respectively. 
These  are  also  shown  in  Fig.  50.  In  the  westerly  latitudes  these  latter 
values  are  much  larger  than  in  the  tropics.  It  will  be  assumed  that 
these  values  are  latitudinally  representative  of  the  oceans  even  though 
the  data  does  have  longitudinal  bias— being  concentrated  along  the  ship¬ 
ping  lanes.  This  is  not  felt  to  be  a  very  restrictive  assumption.  Figs. 
51  and  52  portray  the  vertical  variation  of  wind  speed  and  kinetic  energy 

up  to  2  km  height  by  latitude. 

Stress  Determinations.  A  number  of  researchers  have  attempted 
direct  determinations  of  the  bulk  aerodynamic  drag  coefficient  (C^) 
over  the  oceans  in  various  locations.  Estimates  range  from  1  to 
2  x  10~3.  Brocks  and  Krflgermeyer  (1970)  and  Hasse  (1970)  have  re¬ 
cently  discussed  and  presented  new  data  and  interpretation  on  the  drag 
coefficients  under  neutral  conditions.  They  obtain  a  value  of  igi  1.  3  x 

O 

10  .  Fig.  53  (from  Brocks  and  KrOgermeyer)  shows  no  significant 

change  of  Cj  with  wind  speed.  It  will  thus  be  assumed  that  Cj  is  con¬ 
stant  with  both  wind  speed  and  latitude,  even  though  other  evidence  is 
available  to  indicate  a  slight  increase  of. with  wind  speed. 


difference  is 
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-3 

Given  this  value  of  (ie.  1.3x10  )  and  the  observed  average  of 

the  square  of  the  wind  speeds  at  ship  deck  level  (uQ  ),  the  surface  stress 
<TXZc  ()  can  be  obtained  from  equation  (1)  under  the  assumption  of  neu¬ 
tral  lapse  conditions.  Thus,  with  PQ  -  surface  density. 


C  .p  u 

d  o  o 


(1) 


The  latitudinal  distribution  of  surface  stress  calculated  this  way  is  por¬ 
trayed  by  the  top  curve  of  Fig.  54.  Values  range  from  0.  7  to  2 
dyn/cm^.  These  estimates  of  stress  agree  quite  well  with  the  previous 
estimates  of  Priestley  (1951),  Hellerman  (1967),  Hantel  (1970),  and 
the  estimates  of  others  as  reported  by  Malkus  (1960),  and  Roll  (1965). 
The  bottom  curve  of  this  figure  portrays  the  stress  calculationed  by  the 
geostrophic  departure  method  defined  by  equation  (2) 


lkm 

pi  (v~vg) 


(2) 


'  sf  c 

where  f  is  the  Coriolis  parameter, 

(v-Vg)  is  the  ageostrophic  or  perpendicular  component  of  the 

wind  (  positive  to  the  left  of  the  wind  looking  downstream), 
and 

Si  is  assumed  to  be  1  km  (ie.  no  ageostrophic  wind  at  1  km). 
These  latter  stress  values  are  but  10  to  50  percent  of  the  stress  calcu¬ 
lated  from  (1).  Equation  (2)  grossly  underestimates  the  expected  oceanic 
surface  stress,  especially  in  tropical  latitudes.  The  observed  ageostro¬ 
phic  component,  due  to  frictional  veering  of  only  10°,  is  too  small  to 
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account  for  the  expected  values  of  stress.  Over  land,  where  the  typ¬ 
ical  frictional  veering  is  w  25-35°,  equation  (2)  gives  a  much  more 
realistic  stress  approximation. 1 

Surface  Ageostrophic  Wind  Components.  Fig.  55  shows  the  latitude 
distribution  of  the  surface  or  10  meter  ship  deck  wind  (u)  and  the  cal¬ 
culated  geostrophic  wind  (Ug).  The  ratio  of  u/ug  ranges  from  .  76  to 
.90.  These  are  comparable  to  the  estimates  of  other  researchers.  In 
arriving  at  the  surface  value  of  Ug  ,  a  constant  thermal  wind  was  as¬ 
sumed  in  the  first  and  second  km  layers  for  westerly  winds  and  in  the 
first  and  second  500  meter  layers  in  the  trade  wind  belt.  Ug  is  obtain¬ 
ed  by  the  downward  extrapolation  of  wind  velocity  from  higher  levels 
under  assumed  constant  baroclinicity.  This  will  be  a  good  approxima¬ 
tion  if  no  significant  1st  and  2nd  km  (or  0-0.  5  and  0.  5-1.0  km  for  the 
trades)  thermal  wind  differences  are  present.  In  the  statistical  aver¬ 
age  this  is  felt  to  be  a  reasonable  assumption. 

Fig.  56  portrays  the  latitudinal  distribution  of  the  difference  be¬ 
tween  the  actual  and  geostrophic  surface  wind  components  along  -(u-u_) 

O 

and  perpendicular  (v-Vg)  to  the  boundary  layer  pressure  gradient.  These 
ageostrophic  components  have  similar  magnitudes  and  change  little 
with  latitude. 
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Land  vs.  ocean  veering  is  discussed  in  a  later  section. 
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AGEOSTROPHIC  PARALLEL  (u)  AND 
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Fig.  56.  LATITUDE 

Latitude  distribution  of  the  comparison  of  the  ageostrophic  parallel  -(u-ug) 
and  perpendicular  (v-v^)  wind  component  measured  at  the  10  meter  ship  height. 


Eq.  10*  20*  30*  40*  50*  60*  P 

Fig.  57.  LATITUDE 


The  mean  lowest  km  accelerations  along  the  direction  of  wind  motion  if  zero 
stress  is  present  at  1  km  level.  Note  that  t lie  friction  or  deceleration  term  is 
much  larger  than  the  acceleration  term.  Most  of  this  difference  ^-is  believed 
to  be  accounted  for  by  momentum  convergence  from  upper  layers. 
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MODEL  OF  ZERO  TURBULENT  STRESS  AT  1  km  (  ^  900  mb) 


In  this  section  it  is  assumed  that  the  turbulence  stress  and  ageo- 
strophic  wind  component  decreases  to  zero  at  1  km.  These  assump¬ 
tions  are  made  for  simplicity.  They  lead  to  an  underestimate  of  the 
agoostrophic  wind  component  and  of  the  kinetic  energy  generation  if 
the  top  of  the  boundary  layer  is  at  a  higher  level. 

Boundary  Layer  Acceleration  for  Zero  Stress  at  1  km.  With  the 
observed  values  of  veering  and  stress  obtained  from  equation  (1),  it  is 
possible  to  calculate  the  total  time  derivative  of  the  wind.  Assuming 
no  net  curvature  influences,  the  tangential  equation  of  motion  can  be 
expressed  as 


du 

[mean  through 

1 

35  - 

f (v— v  )  6z  +  -~ 

dt 

[  1st  km 

Az 

g  Az 

_1  xz 
p  3z 


6z 


(3) 


Fig.  57  shows  the  mean  1st  km  frictional  deceleration 
along  the  direction  of  flow  and  the  observed  boundary  layer  accelera¬ 
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from  cross -contour  flow  computed  from  the 
observed  values  of  uQ,  C^,  v,  v^  and  with  A  z  equal  to  assumed  bound¬ 
ary  height  of  1  km.  A  large  (*-*»  5  m/sec  per  day)  deceleration  (-  ^-) 
is  obtained.  This  is  larger  than  expected.  For  steady  conditions  the 
importance  of  a  higher  level  momentum  source  and  downward  transport 
mechanism  is  clearly  evident. 

Kinetic  Energy  Dissipation-Generation  for  Zero  Stress  at  1  km. 


Boundary  layer  kinetic  energy  (KE)  dissipation  per  unit  mass  can  be 
approximated  by. 


1  km 


(4) 


1  f  'az 

KE  Dissipatioi  =  ^  J  U  •  ^7  0  «* 

o 

where  is  taken  to  be  1  km  and 

C  .  is  1.3  x  10"3 
d 

u  is  the  wind  along  the  direction  of  flow. 

Assuming  that  the  stress  decreases  with  height  in  proportion  to  the 
wind  veering,  the  actual  vertical  distribution  of  dissipation  can  be  ob¬ 
tained. 

1st  km  kinetic  energy  (KE)  generation  per  unit  mass  may  be  esti¬ 
mated  from  the  formula 

1  km 

KE  generation  fu  (v-vg)  6Z  (5) 

sfc 

Latitudinal  distributions  of  lowest  km  net  kinetic  energy  (KE)  dissipa¬ 
tion  and  generation  for  all  wind  .  lasses  are  given  in  Fig.  58.  Val- 
ues  are  portrayed  in  watts/m  .  Note  the  large  net  oceanic  dissipation 
to  generation  ratios,  especially  in  the  Tropics.  This  much  larger  than 
one  (dissipation  to  generation)  ratio  is  valid  for  all  speed  classes.  Figs. 
59-62  show  this  same  ratio  for  wind  speed  stratifications  (based  on  the 
500  m  wind)  in  the  ranges  of  0-4,  5-7,  8-11,  >12  m/sec.  For  all 
wind  s-  wUS,  dissipation  is  much  larger  than  generation.  This  is  dif¬ 
ferent  then  over  land  areas  where,  as  reported  by  Kung  (1967)  in  a  five 
year  study  over  North  America,  the  KE  dissipation  and  generation  are 
close  to  balancing. 
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KINETIC  ENERGY  DISSIPATION- GENERATION 
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Fie.  58.  Latitude  distribution  of  the  kinetic  energy  dissipation 
and  generation  for  all  wind  speed  values  if  zero  stress 
is  assumed  at  1  km  height. 
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Fig.  59.  Same  as 


Fig.  58  but  for  0-4  m/sec  winds. 
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LATITUDE 

Fig.  60.  Same  as  Fig.  58  but  for  5-7  m/sec  winds. 
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Fig.  61.  Same  as  Fig.  58  but  for  8-11  m/sec  winds. 
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LATITUDE 

Fig.  62.  Same  as  Fig.  58  but  for  wind  >12  m/sec. 


The  oceanic  dissipation  to  generation  ratio  is  especially  large  in 
the  tropical  latitudes  where,  (due  to  small  value  of  f  ),  pressure  gra¬ 
dients  are  small.  A  typical  8-10°  surface  cross  contour  flow  in  the 
tropics  will  generate  much  less  KE  than  a  similar  crossing  angle  at 
higher  latitudes.  Table  1 1  portrays  idealized  latitudinal  KE  dissipatio 
and  generation  for  wind  velocities  of  5,  10,  and  15  m/sec.  In  these 
calculations  it  has  been  assumed  that 

1)  frictional  boundary  layer  is  one  km  in  depth  and  the  frictional 
veering  at  the  surface  is  10° 

2)  average  ageostrophic  component  in  boundary  layer^v-Vg^is  one 

half  the  surface  ageostrophic  component,  thus  v-v^  =  1  / 2 
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Table  11 

2 

Theoretical  kinetic  energy  dissipation,  and  generation  (in  watts/m  ) 
ratios  for  different  wind  speeds  from  equations  (4)  and  (5)  for  assumed 
conditions  of  10°  surface  frictional  veering  angle,  =  1.3  x  10~3, 
and  mean  veering  through  lowest  km  (v)  equal  to  one-half  surface  veering. 


LATITUDE 


Eg- 10° 

10-20°  20-30° 

30-40° 

40-50° 

50-60° 

1  v 

1  ON 

|o 

o 

For 

u  =  5  m/sec,  ^\T 

= 

.38  m/sec 

Dissipation  .16 

.16  .16 

.16 

.16 

.16 

.16 

Generation 

.03 

.07  .12 

.16 

,20 

.23 

.26 

Diss./Gen. 

5.30 

2.30  1.30 

1.00 

0.80 

0.70 

0.60 

For  u  =  10 

m/sec,  v  - 

v  *  .85  m/sec 
g 

Dissipation  1.30  1.30 

1.30 

1.30 

1.30 

1.30 

1.30 

Generation  .13  .31 

.51 

.69 

.85 

.99 

1.11 

Diss./Gen.  10.00  4.20 

2.50 

1.90 

1.50 

1.30 

1.20 

For  u  *  15 

m/sec,  5T 

~  = 

1.28  m/sec 

Dissipation  3.37  3.37 

3.37 

3.37 

3.37 

3.37 

3.37 

Generation  .29  .71 

1.15 

1.55 

1.93 

2.25 

2.50 

Diss./Gen.  11.60  4.70 

2.90 

2.20 

1.80 

1.50 

1.40 
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It  is  seen  that  the  dissipation  to  generation  ratio  goes  up  with  wind  speed 
and  is  especially  large  in  tropical  latitudes.  This  is  to  be  expected  where 
dissipation  is  proportional  to  the  cube  of  the  wind  speed  and  generation 
to  fvu.  If  the  frictional  veering  angle  is  approximately  constant  with 
latitude,  then  the  dissipation-generation  ratio  of  u  /fv  must  be  quite 
substantial  at  low  latitudes  where  f  is  small. 

In  a  physical  sense  air  in  motion  in  the  boundary  layer  near  the  equa¬ 
tor  should  be  retarded  as  much  (and  transfer  as  much  momentum  and 
kinetic  energy  to  the  surface)  as  boundary  layer  air  at  higher  latitudes. 
The  smaller  tropical  pressure  gradients,  however,  prevent  a  compen¬ 
sating  generation  as  large  as  that  possible  at  higher  oceanic  latitudes. 
Over  land,  where  the  crossing  angles  are  closer  to  30°  rather  than  10°, 
a  much  closer  (at  least  in  middle -latitudes)  dissipation  to  generation  bal¬ 
ance  is  possible.  Fig.  63  portrays  the  net  1st  km  oceanic  KE  dissipation 
with  latitude  and  a  comparison  with  Kung's  (1967)  five-year  average  KE 
dissipation  over  North  America  at  00Z  and  12Z.  In  middle  latitudes  the 
land  and  oceanic  dissipation  energies  are  comparable. 

Model  of  Zero  Turbulent  Stress  at  2  km  (  v-'~  800  mb).  About  20 
percent  of  the  wind  veering  in  the  lowest  2  km  occurs  in  the  2nd  km 
layer.  If  it  is  assumed  that  the  stress  decreases  as  the  winci  veering 
and  that  the  level  of  zero  mechanical  turbulent  stress  is  at  two  km 
( v-'*  800  mb)  rather  than  one  km  ( 900  mb),  then  the  lowest  km  layer 
dissipation  has  been  overestimated  by  about  20  percent  if  the  decrease 
of  stress  is  directly  related  to  the  magnitude  of  wind  veering.  The  mag- 
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Fig.  63.  Net  kinetic  energy  dissipation  in  the  lowest  km  over  the 
oceans  by  latitude  and  compariosn  with  values  of  Kung  at 
00Z  and  12Z  over  the  North  American  continent. 

nitude  of  the  generation  has  similarly  been  underestimated  by  about 
20%-25%.  This  model  of  the  zero  mechanical  stress  at  2  km  ratherthan 
1  km  leads  to  a  decrease  of  the  deceleration  to  acceleration  and  dis¬ 
sipation  to  generation  ratios  of  about  45%-50%.  Dissipation  to  generation 
ratios  (similar  to  the  previous  figures)  based  on  a  2  km  layer  top  of  the 
boundary  layer  for  the  various  wind  speeds  are  shown  in  Figs.  64-68. 
Even  though  these  ratios  are  decreased,  they  remain  substantially  larger 
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Fig.  64. 
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Fig.  65 


KINETIC  ENERGY  DISSIPATION- GENERATION 
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Fig.  66. 
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KiNETIC  ENERGY  DISSIPATION-GENERATION 
IN  SFC-I.O  Km.  LAYER  IF  ZERO  STRESS 
AT  2.0  Km.  (  >  12  m/sec  WINDS) 
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Fig.  68. 
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Fig.  69.  Kinematic  eddy  viscosity  coefficient  required  for 
observed  stress  gradient  and  vertical  wind  shear 
gradient  in  the  lowest  km  layer. 
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than  one  in  the  topics.  In  the  middle  latitudes,  a  more  balanced  dis¬ 
sipation  to  generation  ratio  is  obtained. 

Ocean  vs.  Land  Frictional  Veering.  Mid-latitude  and  subtropical 
land  areas  show  typical  wind  veering  in  the  1st  km  layer --not  of  8°-10° 
as  observed  over  the  oceans,  but,  as  shown  by  Bonner  and  Smith  (1967), 
Kung  (1967)fand  Gray  (1968)  values  of  25o-30°--about  three  times  larger 
than  oceanic  values.  The  Bonner  and  Smith  study  consisted  of  22  Asian 
and  15  U.  S.  stations  for  winter  and  summer.  12,896  soundings  were 
involved.  The  average  veering  in  the  1  km  layer  was  29.  7°.  The  bound¬ 
ary  layer  kinetic  energy  generation  observations  of  Kung  (1967)  require 
an  average  lower  km  wind  veering  of  25-30°.  The  author's  previously 
undocumented  land  veering  study  consisted  of  a  12  month  composite  of 
the  average  difference  of  crossing  angle  of  winds  with  pressure -height 
lines  at  the  surface  and  850  mb  over  the  eastern  half  of  the  U.  S.  Ap¬ 
proximately  12,000  values  were  involved.  The  average  veering  dif¬ 
ference  between  the  surface  and  850  mb  levels  amounted  to  25-30°. 
See  Table  12. 

Middle  latitude  dissipation  to  generation  ratios  over  land  and  ocean 
are  y  and  ^  respectively.  Are  these  large  land  frictional  veerings 
due  primarily  to  higher  land  drag  coefficient?  Additional  p.  b.  1.  ob¬ 
servational  research  is  presently  going  on  at  Colorado  State  University 
(Hoxit,  1972)  to  try  to  better  understand  these  land-ocean  differences. 

Estimated  Viscosity  Coefficient.  From  the  calculated  stress  of  e- 
quation  (1)  and  the  observed  derivative  of  vertical  shear  in  the  bound- 
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Table  12 

FRICTIONAL  WIND  VEERING  OVER  LAND 


Bonner  and  Smith  (1967) 

22  Aslan  and  15  U.  S.  stations.  Winter  and  summer 
observations,  12,896  soundings. 

Average  veering  in  1  km  layer  29.7° 

Rung  (1967) 

5  years  observations  over  N.  Am.  To  obtain  observed 
kinetic  energy  generation,  frictional  veering  must 
be  in  the  range  of  25°-35°. 

Gray  (1969) 

12-month  composite  of  the  average  difference  of  crossing 
angle  of  winds  with  pressure-height  lines  at  the  surface 
and  850  mb  ht.  Approximately  12,000  observations. 
Average  difference  came  to  v/»  30° . 
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ary  layer,  it  is  possible  (using  500  m  shear  gradients)  to  make  an  es¬ 
timate  of  the  mean  oceanic  p.  b.  1.  kinematic  eddy  viscosity  coefficient 


<r> 


(6) 


Fig.  69  shows  the  latitudinal  estimate  of  this  coefficient  in  units  of 
5  2 

10  cm  /sec.  These  values  are  within  the  central  range  of  previous 


estimates. 
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VII.  DISCUSSION 

Conceptual  View  of  Oceanic  Planetary  Boundary  Layer.  The  author 
views  the  p.  b.  1.  as  being  primarily  composed  of  mechanically-driven 
gust-scale  eddies  ( 50-500m  size)  whose  density  and  magnitude  de¬ 
crease  throughout  the  layer.  Where  cumulus  clouds  are  present,  larger 
scale  eddies  (1-10  km)  can  exist  and  extend  to  high  levels,  as  portrayed 
in  Fig.  70.  The  purely  mechanical  influences  of  the  oceanic  p.  b.  1.  can 
be  viewed  as  typically  extending  up  to  about  one  km  height  ( 900  mb). 


Fig.  70.  Idealized  view  of  the  decrease  of  turbulent  eddies 
with  height  and  without  cumulus  clouds. 
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A  diurnal  p.  b.l.  height  variation  as  is  present  over  land  (due  to  diurnal 
buoyancy  and  wind  variations)  is  not  found  over  the  ocean.  Where  cu¬ 
mulus  convection  (i.  e. ,  moist  processes  are  allowed  for)  oceanic  p.  b.  1. 
can  extend  to  higher  levels.  Cumulus  clouds  can  distribute  surface  mo¬ 
mentum  to  higher  levels  than  the  dry  (and  thus  stable)  purely  mechanical 
turbulence  processes  are  able  to  do.  This  does  not  negate  the  ideal¬ 
ized  view  of  a  mechanical  boundary  layer  of  about  1  km  height,  how¬ 
ever.  Cumulus  clouds  occupy  only  a  small  percent  of  the  ocean  area. 

The  cumulus  clouds'  net  influence  on  altering  the  mechanical  stress 
level  over  the  oceans  (except  when  intense  convection  and  large  .  crtical 
shears  are  present)  is  probably  not  very  great.  For  these  reasons, 
the  concept  of  a  general  one  km  surface  turbulence  boundary  layer  is 
felt  to  be  "statistically"  valid. 

Required  Downward  Transport  from  Higher  Levels.  It  is  clearly 
seen  (if  these  observations  are  correct)  that  the  oceanic  p.b.  1.  (espec¬ 
ially  in  the  tropics)  is  not  a  self-contained  momentum  and  kinetic  en¬ 
ergy  system.  The  lowest  km  layer  is  continually  running  itself  down. 
Replacement  kinetic  et  ergy  must  come  from  higher  levels,  not  by  gust- 
scale  transfers  but  by  meso  and  synoptic  scale  clear  air  sinking  motions, 
as  pictorially  visualized  in  Fig.  71,  or  by  some  other  processes  such 
as  cumulus  (or  Cb)  downdraft  (Zipser,  1969)  transfers.  It  may  also  be 
possible  for  vertical  transports  to  be  accomplished  by  cumulus  (or  Cb) 
induce^.  gravity  waves. 


cumuiut-tcol* 
boundary  lay «r 


VERTICAL  TRANSPORT  Q£  HORIZONTAL  MOMENTUM 
IQ  OCEAN  SURFACE. 


By  gust  scole 


By  cumulus-  meso-  and 
synoptic  scales 


By  gust  and 
cumulus  scales 


I'ig.  71.  Idealized  view  of  the  downward  transfer  of  mom¬ 
entum  by  the  different  scales  of  motion. 
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Fig.  72.  Required  latitudinal  average  up-moist  and  down-drv 
vertical  circulation  at  the  1  km  level  required  to  ex¬ 
plain  the  observed  rainfall  and  the  net  heat  losses  by 
radiation  and  evaporation  as  discussed  by  Gray  (1 ! ‘ 7 1> . 
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These  larger  scale  downward  vertical  transfers  are  believed  to  oc¬ 
cur  in  association  with  the  required  up-moist  and  down-dry  vertical  cir¬ 
culations  necessary  to  produce  the  rainfall  and  balance  the  net  lower 
tropospheric  radiational  and  evaporational  cooling.  The  author  (Gray 
1971)  has  previously  discussed  this  required  up-moist  and  down-dry 
vertical  circulation  whose  magnitude  in  mb/day  is  shown  in  Fig.  72. 

The  required  down-dry  or  clear  air  sinking  will  replace  a  substantial 
portion  of  the  p.  b.  1.  air  every  day,  especially  in  tropical  and  subtrop¬ 
ical  latitudes.  Fig.  73  is  an  idealized  view  of  the  typical  kinetic  energy 
balance  occurring  at  tropical  and  middle  latitude  locations.  This  was 
determined  as  a  compromise  between  assuming  a  one  and  two  km  zero 
stress  height. 

Depth  of  the  Boundary  Layer.  The  self-contained  or  steady  bound¬ 
ary  layer  theories  (Ekman,  1905  or  Charney,  1969)  require  that  the 
p.  b.l.  frictional  dissipation  be  balanced  by  an  erual  cross  contour  flow. 
Where  pressure  gradient  becomes  small,  as  where  one  approaches  the 
equator,  these  theories  require  the  depth  of  the  l  oundary  layer  to  in¬ 
crease  so  the  weaker  cross-contour  flow  can  ac  t  over  a  larger  vertical 
depth  and  still  balance  the  surface  dissipation.  The  observations  sum- 
marized  in  this  report  do  not  show  a  general  increase  in  the  thickness 
of  the  boundary  layer  as  one  approaches  the  equator.  Rotation  appears 
not  to  be  a  crucial  parameter  in  determining  the  depth  of  the  p.  b.  1. 

And  why  should  it?  Tne  author  believes  that  the  depth  of  the  boundary 
layer  is  primarily  determined  by  the  intensity  of  gust-scale  kinetic 
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energy  and  by  the  vertical  stability.  These  factors  are  not  related  to 
latitude. 

If  the  p.  b.  1.  is  not  a  self-contained  momentum-energy  system,  there 
is  no  reason  why  a  systematic  latitude  variation  need  occur.  Why  should 
the  mode  of  turbulent  interface  momentum  transfer  from  atmosphere 
to  surface  vary  with  latitudes?  Do  the  turbulence  eddies  know  at  what 
latitude  they  exist? 

Vertical  lapse  rates  in  the  ocean  planetary  boundary  layer  are  typ¬ 
ically  0.6-0.  8  of  the  dry  adiabatic  values.  This  means  that  the  turbu¬ 
lence  eddies  (rising  at  the  dry  adiabatic  rate)  must  continually  act  a- 
gainst  a  stable  environment.  The  higher  they  go,  the  more  negative  is 
their  buoyancy.  This  vertical  stabilizing  acceleration  of  the  environ¬ 
ment  (Az)  may  be  expressed  as  a  function  of  height  (z)  as 


A 

z 


sf  c 


(yd  -  ya) 
T 


6z 


(7) 


where 


yd,  ya  = 


z 


T 


dry  and  actual  lapse  rates 

top  of  p.  b.  1.  (level  at  which  u'w'  approaches 
zero,  where  u'  and  w'  are  horizontal  and  vertical 


turbulent  eddies. 
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Assuming  the  following  mean  layer  values  of  actual  lapse  rate  (Ya  ) 
and  average  temperature  (T ) 


Layer  (km) 

<Ya) 

T 

0  -  0.5 

•8Yd 

290 

0.5  -  1.0 

•  7Yd 

287 

1.0  -  1.5 

•  8Yd 

284 

1.5  -  2.0 

•  6  Yd 

281 

the  mean  value  of  (in  cm/sec  per  gm)  through  these  layers  is 


Layer 

A 

z 

0-0.5 

km 

1.8 

0. 5-1.0 

km 

6.4 

1.0-1. 5 

km 

13.0 

1. 5-2.0 

km 

20.6 

The  vertical  stabilizing  energy  may  then  be  expressed  as 


E 

z 


(vertical  stabilizing  energy) 


■  Z>j> 

Az  6z 
sf  c 


(Z_  -  Z  ) 

T  src 


(8) 
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c  9  2 

The  values  of  E,  (in  10  cm  /sec  per  gm)  at  various  heights  (for 
z 

these  Az  values)  are 

Height  (km) 


0.  5 

0.9 

1.0 

4.  1 

1.5 

10.6 

2.0 

20.9 

It  is  often  observed  that  gustiness  of  surface  wind  speeds  causes 

fluctuations  from  about  half  to  one-and-a-half  times  the  mean  wind 

speed.  The  typical  fluctuation  of  wind  with  mean  speed  of  5,  10,  and 

15  m/sec  may  be  thought  of  as  being  in  the  range  of  about  2-1/2  to 

7-1/2  m/sec,  5  to  15  m/sec,  and  7-1/2  to  22-1/2  m/sec  respectively. 

Assuming  (for  rough  magnitude  considerations)  that  the  typical  gust- 

scale  horizontal  turbulent  kinetic  energy  (l/2  u'  ,  where  u'  is  a  wind 

eddy  and  the  bar  —  indicates  a  long  time  average),  is  equal  to  one- 

,  "  2 

half  and  one  times  the  kinetic  energy  (KE)  of  the  mean  flow  or  (l/2u' 

9  2  5  2  2 

1/4  u  and  l/2  u  ),  then  the  turbulent  kinetic  energy  (in  10  cm  /sec 

per  gm)  is  for  various  mean  wind  speeds: 

Turbulent  KE 


u  (m/sec) 

if  1  /4  “  “ 

if  1/2  u1 

5 

CO 

• 

o 

1.2 

10 

2.  5 

3.  0 

15 

5.  6 

11.2 

20 

10.0 

20.  0 

-  VH  ^  •  - 
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Assuming  the  mixing  length  hypothesis  to  be  valid  to  the  extent  that 
the  horizontal  (1/2  u'  )  and  vertical  (  1/2  w'  )  kinetic  energies  are 
about  the  same,  one  might  equate  the  horizontal  turbulent  kinetic  energy 
to  an  equal  restrictive  buoyant  energy  height;  i.  e„,  E ^  (ht  /  trrbu- 
lent  KE.  This  height  will  be  called  the  top  of  the  mechanical  turbulent 
p.  b.  1.  It  is  the  height  above  which  surface  generated  wind  eddies  can¬ 
not  penetrate  due  to  buoyant  damping.  When  this  is  done  for  the  above 
assumed  kinetic  energies,  the  following  stabilizing  restrictive  heights 
are  obtained: 


u  (m/sec) 


o 

10 

15 


Height  (km) 
Turb.  KE=l/4"u2 


0.4 
0.7 
1.  1 


Turb.  KE  —  i/2  u" 2 


0.  7 
1.  1 
1.5 


In  active  convective  situations  where  the  surface  wind  speeds  are 
strurg  and  increase  sharply  through  the  troposphere,  the  top  of  the  p.  b.  1. 
should  be  expected  to  go  higher  than  1  km.  When  wind  speeds  decrease 
with  height  in  the  lower  troposphere,  as  in  the  trade  winds,  the  u'w' 
correlation  in  cumulus  clouds  will  be  of  opposite  sign  to  the  surface 
u'w'  eddy  correlation,  and  a  cross-over  or  zero  stress  height  should 
be  readily  observed  at  low  levels. 

Conclusion- -The  height  of  the  p.  b.  1.  can  be  specified  in  a  general 
sense  by  vertical  stability  and  turbulence  intensity  arguments  by  them¬ 
selves  without  resort  to  rotation  or  steady  state  assumption.  The  in- 
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formation  of  Figs.  40-43  does  indicate  an  increase  of  the  mechanical 
boundary  layer  thickness  with  increase  of  wind  speed.  Other  data  cur¬ 
rently  being  evaluated  by  Hoxit  (1972)  over  land  also  are  showing  sim¬ 
ilar  p.  b.  1.  height  variations  with  similar  stability  and  wind  speed 
changes. 

Importance  of  Oceanic  Boundary- Layer  Frictional  Veering  and  Re¬ 
lation  to  CISK  Mechanism.  The  magnitude  of,  and  relevance  of,  the 
first  km  frictional  wind  veering  over  the  oceans  has  been  open  to  some 
misinterpretation.  Charney  and  Eliassen  (1964),  Ooyama  (1964),  and 
the  author  Gray  (1968)  have  proposed  it  as  part  of  an  important  phy¬ 
sical  mechanism  for  the  intensification  of  cloud  clusters  and  tropical 
storms.  Charney  and  Eliassen  (1964)  have  chosen  to  call  the  low-lev¬ 
el  frictional  forcing  and  the  meso-scale  flow  instability  it  produces 
(applied  to  tropical  storm  development)  Conditional  Instability  of  the 
Second  Kind  (CISK).  The  CISK  idea  has  been  used  as  the  physical  basis 
for  explaining  tropical  storm  development  and  has  been  extensively  em¬ 
ployed  in  hurricane  intensification  models. 

Some  meteor  ologists,  on  the  other  hand,  who  believe  in  the  primary 
importance  of  baroclinic  as  opposed  to  frictional  forcing  processes  have 
questioned  the  dominant  role  of  the  CISK  instability  idea.  They  feel 
that  it  has  been  overly  stressed.  This  paper  does  not  concern  itself 
with  the  entire  CISK  proces  s  but  only  with  that  part  dealing  with  low- 
level  frictional  forcing. 
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The  CISK  idea  implies  the  dual  physical  processes  of  frictional  con¬ 
vergence  and  unstable  feedback  growth.  But  this  kind  of  cumulus -broad¬ 
er  scale  cooperative  instability  only  takes  place  with  intensifying  trop¬ 
ical  clusters  and  storms  —not  a  frequent  phenomenon.  The  usual  oc¬ 
currences  of  frictionally-forced  cloudiness  along  frontal  systems,  with 
cloud  clusters,  the  ITCZ,  etc.,  takes  place  without  broader-scale  flow 
intensification.  The  CISK  instability  connotation  is  not  an  accurate 
description  for  the  frictionally-forced  convergence  by  itself.  For  sim¬ 
plicity  we  might  call  this  low-level  frictionally-forced  process  by  itself 
"Low-level  Inflow  from  Frictional  Turning  (LIFT).  This  low  level  fric¬ 
tional  LIFT  is  only  part  of  CISK,  thus 

CISK  =  LIFT  +  (Feedback  Instability) 

CISK  implies  both  the  frictional  forcing  and  the  cumulus  feedback  insta¬ 
bility.  It  is  important  to  note  this  difference.  Unique  tropospheric  flow 
features  above  the  p.  b.l.  must  be  present  (i.  e. ,  low  tropospheric  vert¬ 
ical  wind  shears,  conditionally  unstable  lapse  rates,  high  middle  level 
moisture  content,  etc.)  for  cumulus  feedbacks  and  broader-scale  flow 
intensification  to  occur.  Except  in  the  comparatively  rare  cases  of  in¬ 
tensifying  clusters  and  storms,  this  feedback  instability  potential  is  not 
released.  This  should  not  be  taken  as  a  reason  for  denying  the  physical 
importance  of  the  LIFT  process  by  itself,  however. 

Importance  of  Frictionally-forced  Convergence.  The  author  believes 
this  LIFT  mechanism  (at  least  the  physical  idea  of  it)  to  be  a  1'undamen- 
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Fig.  74.  Portrayal  of  how  cyclonic  shear  or  relative  vorticity  in  a 
zonal  non-divergent  current  at  950  mb  can  produce  p.  b.  1. 
convergence  if  a  frictional  veering  of  10°  were  present. 

It  is  in  these  shear  areas  when  most  deep  cumulus  convec¬ 
tion  takes  place. 
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tally  important  one  over  the  oceans  in  specifying  the  location  and  den¬ 
sity  of  cumulus  convection.  A  frictional  wind  veering  in  the  lowest  km 
of  10°  would  lead  to  surface  cross -isobaric  flow  (sin  10°  -  .17)  of 

10%-20%  of  the  total  wind.  The  average  lowest  km  perpendicular 
component  would  be  8%-10%  of  the  wind.  In  an  average  over  the 
oceans,  then,  frictional  forced  lowest  km  convergence  (C^,)  would  occur 
where  there  is  horizontal  wind  shear  as  in  Fig.  74  and  would  be  related 
to  the  relative  vorticity  (Zy  )  as 

Ct  ~  ^  ty  m 

As  most  oceanic  cumulus  clouds  have  their  bases  at  v^l  km  level,  a 

lowest  km  forced  convergence  of  10  percent  of  the  relative  vorticity 
would  be  a  very  significant  cumulus  production  mechanism  in  regions  of 
large  horizontal  wind  shears  such  as  exist  along  frontal  zones  and  on 
the  equatorial  side  of  the  Trades.  It  is  no  mere  coincidence  that  it  is 
in  these  regions  where  the  majority  of  the  deep-cumulus  are  located. 

Conclusion- -There  is  a  significant  (but  not  large)  frictional  wind 
veering  in  the  p. b.  1.  ov  er  ihe  oceans.  An  important  purpose  of  the 
GATE  and  GARP  programs  should  be  a  more  exact  quantitative  specifi¬ 
cation  of  this  physical  process. 

Further  Work.  The  author  will  next  attempt  to  stratify  this  observa¬ 
tional  data  by  longitude  and  vertical  stability  and  will  investigate  veering 
at  higher  levels.  He  will  also  obtain  a  larger  statistical  sample  on  the 
veering  closer  to  the  equator.  Extensive  data  stratifications  are  also 
proceeding  over  land  areas  so  that  better  land-ocean  comparisons  can 


be  made. 
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